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ABSTRACT 


Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is 
determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion 
of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We 
explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild 
black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and 
conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We 
test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using 
ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the 
ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization 
evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen 
parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon 
Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid 


velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. 


Keywords: Accretion (14), Black holes (162), Polarimetry (1278), Magnetic Fields (994) 


1. INTRODUCTION 


The Event Horizon Telescope (EHT) Collaboration has re- 
cently published the first images of a black hole (Event Hori- 
zon Telescope Collaboration et al. 2019a,b,c,d,e,f, 202 1a,b; 
hereafter EHTC 1-УШ). These images achieve a diffraction- 
limited angular resolution that corresponds to approximately 
5С М/с?, where M is the mass of the black hole. They reveal 
a bright ring of emission with a twisting polarization pattern 
and a prominent rotationally symmetric mode. 

The polarization structure in the EHT images depends on 
details of the emitting plasma, principally the magnetic field 
geometry. However, it is also affected by the strongly curved 
spacetime near the black hole. Over the past few decades, 
simulated polarimetric images of black holes have been stud- 
ied as a means to understand astrophysical properties of 
their surrounding accretion flows (e.g., Bromley et al. 2001; 
Shcherbakov et al. 2012; Moscibrodzka et al. 2017; Jiménez- 
Rosales & Dexter 2018; Palumbo et al. 2020) and to infer the 
disk inclination and black hole spin through the effects of par- 
allel transport (e.g., Connors et al. 1980; Broderick & Loeb 
2006; Li et al. 2009; Schnittman & Krolik 2009; Gold et al. 
2017; Marin et al. 2018). 

While they are becoming increasingly realistic, these sim- 
ulations are generally difficult to use for broad parameter sur- 
veys because of their computational cost, and they often pro- 
vide little insight into how to decouple astrophysical and rel- 
ativistic effects. 

In this article, we develop a simple toy model to understand 
polarimetric images of black holes. This model consists of 
a ring of magnetized fluid orbiting a Schwarzschild black 
hole. Our model allows arbitrary emission radius, magnetic 
field geometry, equatorial fluid velocity, and observer incli- 
nation. With a single approximation, described in section 2, 
we can analytically compute the resulting polarimetric image 
and can assess its dependence on the input parameters. 

In section 2, we describe the toy ring model and work out 
the relevant relativistic transformations from the frame of a 
radiating fluid element in the ring to the image as seen on the 
sky by an observer. In section 3, we present a series of ex- 
amples to illustrate the primary model features. In section 4, 


we provide analytic estimates of image diagnostics — the ap- 
parent shape of the ring, the vector polarization, and the co- 
efficient of rotational symmetry (85; Palumbo et al. 2020). In 
section 5, we discuss the suitability of our model for com- 
parisons with observations, focusing on the EHT images of 
М87* and polarization “loops” seen during flares of Sagittar- 
ius A* (Sgr A*). In section 6, we summarize our results. 


2. THE MODEL 


We consider an accretion disk around a Schwarzschild 
black hole of mass M. We use standard geometrized units: 
G = с = 1. The fluid radiates from the equatorial plane 
within a narrow range of radii centered on a dimensionless 
radius А, measured in units of M (or GM/c?, including the 
physical constants). With respect to a distant observer, the 
ring is tilted from а face-on orientation by an angle 0j. We 
assume that the tilt is towards the North, so that the line-of- 
nodes between the ring orbital plane and the observer's sky 
plane is in the East-West direction. We take the sky angu- 
lar coordinate x to be oriented towards the West (1.e., to the 
right), and the coordinate y towards the North (i.e., towards 
the top). The fluid has radial and tangential components of 
velocity in the plane of the ring, but no vertical velocity. In 
the comoving frame of the fluid, the magnetic field has ra- 
dial, azimuthal and vertical components. For simplicity, we 
assume that both the velocity field and the magnetic field are 
axisymmetric, though the equations developed in this section 
are valid even without this assumption. 

We wish to compute the following primary observables: 
(1) the shape of the ring as viewed by the distant observer, (2) 
the variation of the polarized intensity around the observed 
ring, and (3) the orientation and pattern of the polarization 
vectors around the ring. An exact calculation requires inte- 
grating the geodesic equation, which has to be done numer- 
ically. However, with one simplification, described below, it 
is possible to do all the calculations analytically. This sim- 
plified model provides a convenient method for investigating 
polarization properties of idealized models. 


2.1. Geometry, Lensing and Special Relativity 
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Figure 1. Geometry in the geodesic frame, or G-frame. In the 
Schwarzschild metric, each null geodesic is confined to a plane that 
intersects the black hole. The G-frame, defined for photons emitted 
at point P and reaching a distant observer at relative angle y, corre- 
sponds to Cartesian axes centered on the black hole, with $ in the 
direction of P and the $-2 plane given by the geodesic plane. We 
approximate the emission angle a in this frame using Equation 4. 


In the ring plane, we consider a fluid element P located at 
azimuthal angle ф measured from the line-of-nodes. We are 
interested in a null geodesic, a light ray, that travels from P 
to the observer. This geodesic lies in a plane that includes 
the line from the black hole O to the point P, as well as the 
line from O to the observer (see Fig. 1). We set up Cartesian 
coordinates in the geodesic plane so that the unit vector along 
the x-axis ĉis oriented along OP and the observer lies on the 
$-2 plane. We call this the geodesic frame, or G-frame. The 
angle y between X and the unit vector fi towards the observer 
satisfies 


cos у = — sin 6, sin ф, 
sin y = (1 — cos? y)! ^. (1) 


We consider a null geodesic with conserved energy! К, = 
—1 traveling from Р to the observer. At the location P, 
the orthonormal time component kio of its 4-wavevector 
is given by (the redshift factor here is calculated using the 
Schwarzschild metric, as appropriate for the assumed non- 
spinning black hole) 

? К, 1 
Ko" Ea" Gay ©) 


where the subscript “(С)” indicates that this quantity is mea- 
sured in the G-frame. Also, since the geodesic lies in the 


xz-plane, we have E — 0. To determine the other two com- 
ponents of k, we need the angle a in Fig. 1, in terms of which 


we can write 
£ _ yt 2 Hd os 
Ka) = ko) cosa, KG) = kg) sin a. (3) 


! This is the photon energy measured by an observer at infinity, and we 
normalize it to unity. 


Line of Nodes 


Figure 2. Geometry in the P-frame. This frame is aligned with the 
rotating gas at emission radius А and emission azimuth ¢. The £ĉ 
direction lies along the radial line from the black hole at O to the 
emission point P, and $ is the azimuthal direction. The equatorial 
magnetic field Beg and fluid velocity B lie at angles 7 and y to € 
in the x-y plane, respectively. Our model allows these angles to be 
specified independently, but we will later focus on the physically 
motivated choices of 7 = y and ņ = y + л (see section 3). 


Instead of attempting to calculate o precisely, which would 
require a numerical integration of the geodesic equation, 
we use the following approximate formula obtained by Be- 
loborodov (2002), 


2 
COS œ = cosy + в“ — cos y), 


sina = (1 - cos? a)! ?. (4) 


This approximation is surprisingly accurate even for values 
of R of order a few (see sec. 3.6 and Appendix A). 

We now switch to a Cartesian frame that is aligned with the 
orbiting fluid ring. We take € along ОР, ӯ in the azimuthal di- 
rection at P parallel to à, and 2 perpendicular to the orbital 
plane. We call this the P-frame (see Fig. 2). The G-frame 
and P-frame have a common X-axis. Therefore, transform- 
ing from one to the other involves rotation by some angle £ 
around the x-axis. To determine £, we note that the unit vec- 
tor f; from the black hole О towards the observer has Carte- 
sian components (cos y, 0, sin y) in the G-frame, and Carte- 
sian components (— sin Ө, sin ф, —sin 0, cos ф, cos 6») in the 
P-frame. Since a rotation by angle £ transforms one set of 
components to the other, we obtain 


0, ; in o 
eee cos , dap sin Ө, cos A (5) 
sin y sin y 
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Applying this rotation to the orthonormal components of Кс), 
we obtain the corresponding orthonormal components in the 
P-frame, 


» 1 М cosa 


5 sin é sinc 5 cos & sina 
Ф) = Gay (7 Ea О) 


The fluid at the point Р moves in the xy-plane of the lo- 


cal P-frame with a velocity B. which we write in the local 
Cartesian coordinate frame as (see Fig. 2) 


B = B (cosy $ + siny 5) . (8) 


Our sign convention is that radial motion towards the black 
hole corresponds to cosy < 0, and clockwise rotation on 
the sky corresponds to siny < 0. In the case of M87*, the 
rotation is clockwise. The velocity B describes motion of the 
fluid through the ring; the ring model itself is not expanding 
or contracting. 

We now transform to the fluid frame — the F-frame — 
by applying a Lorentz boost with velocity B. This gives the 
following orthonormal components of k, 


kip = y kp) — УВ cosx kee) — yB siny kp 


Кё = —yB cos y Kip + (1 + (y — 1) cos? y) К}, 
+(y — 1) cos y siny kp 

Kp = УВ sin y kip + (у — 1) cos y siny Кф, 
+(1 + (у = Dsin? y) is 


Ky = Юр): (9) 


2.2. Transformation of Polarized Intensity 


Any radiation emitted along ks in the F-frame is Doppler- 


shifted by the time it reaches the observer. Since Ко in the 
observer frame is equal to unity, the Doppler factor 6 is 
Ko 1 
ô = Га = "ND 
(F) (F) 

This includes both gravitational redshift and Doppler shift 
from velocity. 

In the fluid frame, there is a magnetic field which we write 
аз“ 


(10) 


B=B,3+ BS + BZ 
= Beg (cos  $ + sing $) + В, 2 (11) 


z A 
= Beg + В, 


2 Because the emission of synchrotron radiation is best described in the 
fluid frame, we find it convenient to specify the magnetic field components in 
this frame. The %, $, 2 axes in the fluid frame are related to the corresponding 
axes in the P-frame (equivalently, the Schwarzschild frame, e.g., eq 19), 
via a Lorentz transformation with velocity B. The transformation of field 
components between the two frames is worked out in Appendix B. 


where the second line describes the field components in the 
equatorial plane in terms of a magnitude Beq and an orien- 
tation 77 (see Fig. 2). The intensity of synchrotron radiation 


emitted along the 3-vector Ер depends on sinZ, where ¢ is 
the angle between Ер апа the magnetic field B: 


> > 
_ [kay x B| 


sin е 
kel |B| 


(12) 


In the case of thermal synchrotron emission, the intensity 
also depends on the ratio of the emitted photon energy Лу to 
the electron temperature kT,. At low frequencies Av « KT,, 
the intensity is proportional to sin?” С (e.g., Mahadevan et al. 
1996), whereas in the opposite limit hy > КТ», the inten- 
sity varies as a very large positive power of sin Z, because of 
the exponential cutoff of the particle energy distribution and 
the corresponding rapid decline of emissivity with increas- 
ing frequency. In general, if the emitted intensity varies as 
I, ~ v^*5, then the angle dependence goes as (sin Z)!*^». In 
models of M87*, a dependence ~ sin? £ is often obtained at 
230GHz. This corresponds to a, ~ 1, which is consistent 
with the synchrotron emission being close to its peak at this 
frequency (vF, roughly constant). In the analysis below, we 
explicitly retain the a, dependence. However, we set a, = 1 
for the numerical calculations described in sec. 3, and also 
when we series-expand the equations in Appendix D. 

The factor (sin ё)!+% discussed in the previous paragraph 
is the emission per unit volume. To convert this to the emerg- 
ing intensity in the fluid frame we need to multiply by the 
geodesic path length / through the emitting region. We as- 
sume that the medium is optically thin to its own emission. 
If we model the emitting fluid as a thin disk of vertical thick- 
ness H, then the path length is 

кі 


b= ӧн. (13) 


So far, we have discussed the emitted intensity in the fluid 
frame. This intensity is Doppler-boosted by a factor of 6°*% 
by the time it reaches the observer? Thus, the intensity |P| 
of linearly polarized synchrotron radiation that reaches the 
observer from the location P is 


|P] = ét lp [Bite зїп!+® ¢ (14) 
61,182 sin? £ for a, = 1, (15) 


where we have omitted a proportionality constant. Since 18| 
. . . . 9 

is constant around the ring, ће factors involving |B| could be 
eliminated from Equations 14 and 15 and absorbed into the 


3 In the context of a continuous relativistic jet, a Doppler boost factor of 
52+ is generally used (e.g., Blandford & Kónigl 1979). That corresponds 
to the combined quantity / porte, where for motion parallel to the jet axis, 
ly œ 6-1. Our formulation, with / р handled as a separate factor, is more 
general. 
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omitted proportionality constant. We retain these factors be- 

. ә . . . 
cause keeping track of |B|? and its components is convenient 
for much of the analysis in Appendix D.* 


2.3. Transformation of Polarization Vector 
We next work on the polarization vector. In the fluid frame, 
the É-vector of the radiation is oriented along Ер х В, 1.е., 


D 2 . 
perpendicular to both Kp, and B. Therefore, we write the 
orthonormal components of the polarization 4-vector f" as 


; (kæ x B) 
Л 70 fg = Il 
(F) 
(бх). (хд), 
ae, fe аб 
koal kel 


By construction, this 4-vector satisfies 


fk, = 0, f^ f, = sin? СВР. (17) 


An inverse Lorentz boost transforms the 4-vector tes back to 
the P-frame: 


fiy = fi + YB cos x fy + yB sin x frs. 
fie = yBcos y fi + (1 + (y — Dcos? y) fiy 
+(y — 1) cosy siny fr. 
fip = уВзіпд fir) + (y — 1) cos y siny fi, 
+(1 + (у= D sin? y) fi, 
fe = foe (18) 


Since the Cartesian unit vectors £, Ӯ, 2 іп the P-frame аге 
oriented along the spherical polar unit vectors ?, à, — 6, of 
the Schwarzschild frame, the orthonormal components of k 
and f in Schwarzschild coordinates are 


kK = k 


BL 6_ ê à 25 
(p К = Кр k = =k kK =k 


P? e» (9 


F=f flf fe f-f» 00) 
The photon geodesic emitted at P has three conserved quan- 
tities (see for instance Bardeen 1973): its energy k, = —1, 
its angular momentum around the 2 axis ky = КА, and the 
Carter (1968) constant C, which is the square of the total an- 
gular momentum of the photon for the Schwarzschild metric. 
In the P-frame the Carter constant is 


C=R (0) + («) | | (21) 


4 Alternatively, we could assume 18| = 1, as indeed we do in all the 
plots, eliminate IB from Equations 14 and 15, but still keep track of the 
components of Bin Appendix D. 


Using the conservation of kg and C, we find the coordinates 
x and y of the geodesic at the observer sky plane (recall the 
orientation of the sky coordinates x, y described at the top of 
section 2) (Bardeen 1973), 


ke Еф 
sind, — sins’ 


X—— 


cds к)" — cot? 6, ue sen(sing). (22) 


To compute the polarization vector at the observer, 
we make use of the Walker-Penrose constant Кү + iK» 
(Walker & Penrose 1970), which takes a simple form for 
a Schwarzschild spacetime. At the position P, we have (us- 
ing the sign convention in Himwich et al. 2020), 


KQu-R(f'-Kf) К = -ЕЗ( f? – КӘ РФ). (23) 


Both Кү апа К» are conserved along the geodesic. There- 
fore, knowing their values, we can evaluate the two trans- 
verse components of the polarization electric field E at the 
observer. If we use the normalization used in Himwich et al. 
(2020), the field components are 


Exnorm = = 2 , 
nom GI KD 08 € 9)? 
DNE. uc NE 
7 [Ki + К2) О? + ya? 
B un + En =1, (24) 


which is normalized to unity. This normalization is suitable 
for plotting the orientation of polarization vectors in the xy- 
plane. 

An alternative normalization is 


Е = yK + xK; 
yKi — xK5 
Еу= 8 
х= +у 
Е? + Е? = sin? ¢ |B}. (25) 


This retains the original normalization of f^" in the fluid 
frame (eq 17), hence the electric field is proportional to 
sin Z |Ë]. 

For computing the observed polarized intensity, we need 
to include the dependence on the Doppler factor 6 and path 
length /,, and must also ensure the correct powers of sin Z 
and 18| as given in Equations 14 апа 15. Since the intensity is 
proportional to EI ^, we therefore write the observed electric 


8 


field components as 


: E 
E oii = 6090/2 П (sin губ? |8|01+0)/2 Е, 


,norm 
= @Ө+%)/? pa (sin CY DP gie»? Ex, 
(26) 


. E 
Eds = 6(3+2)/2 jun (sin pyrene 18|01+0)/2 Enom 


= g8te2 d (sin е Э?|В®- 7? Е. 
(27) 
Е + Е. E |IP(9), 
where Р(ф) is the observed linear polarized intensity of ra- 
diation that is originally emitted by a fluid element at ring 
azimuthal angle д. 

We need one more transformation: we must convert the 
coordinates (А, 9) of the emitting region in the fluid to the 
Cartesian sky coordinates (x, y), or equivalently the polar sky 
coordinates (p, y), at which the radiation is observed, 


х= рсоѕф, yc-psing. (28) 


The relation between (А, ф) and (р, ф) is worked ош in Ap- 
pendix С. The observed linear polarization Р(ф) can then 
be described in image coordinates by the complex function 
P), 


P(g) = Ay) + UY), (29) 


where the Stokes parameters О(ф) and (е) are obtained 
from the electric field components Е. оьѕ, Ey,ops using Equa- 
tion D.10. The electric vector position angle (or EVPA) is 
then 


1 
EVPA = 5 arctan Е (30) 


This completes the calculation of the intensities О, U, Р 
on the image plane. If one wishes to calculate fluxes in the 
sky plane corresponding to specific source configurations in 
ring coordinates (А, ф), it would be necessary to apply the 
Jacobian of the transformation from (А, ф) to (р, ф), as in 
Figure 10. The Jacobian determinant is evaluated in Ap- 
pendix C. . 

To summarize, in this section we showed how, given 
the position (А, ф, Fig. 2) and velocity (B,y. eq. 8) of a 
synchrotron-emitting fluid element located on a tilted equa- 
torial plane around a Schwarzschild black hole, and given 
also the magnetic field configuration (Вед, 7, Bz, eq. 12) in 
the frame of the fluid, one can calculate the sky coordinates 
(x, y, equivalently р, o) of the image of this radiating ele- 
ment, and the linearly polarized intensity and position angle 
of the observed radiation. The mapping from the radiating el- 
ement to the observer's image plane is written as a sequence 
of analytical calculations that do not require numerically 
integrating the geodesic equation or iteratively solving any 
equation. The equations are written in sufficient detail for 
easy incorporation into modeling calculations. 


3. EXAMPLE MODELS 


The simple model considered in the previous section has 
the following parameters: tilt angle of the ring 6,, ring radius 
R, velocity vector of the fluid B. which is parameterized by 
В = v/c and y (eq 8), fluid frame magnetic field B, which is 
parameterized by either В,, By, Bz, or Beg, 7, В, (eq 12), and 
spectral index оу. Figures 3—5 show the polarization patterns 
produced by this model for selected values of the parameters. 
In all these examples, we choose 0, = 20° and a, = 1. 

Before considering the examples, we briefly summarize 
the salient features of the polarized image of M87* obtained 
by the ЕНТ (EHTC ҮП). First, the linear polarized flux 
shows a pronounced asymmetry around the ring. The po- 
larized flux is strong between PA (measured East of North) 
~ 150° and ~ 300°; the peak polarized intensity is around PA 
200° on April 5 and 240° on April 11. The linear polarized 
flux is much weaker at other angles. The large scale jet in 
М87* is oriented towards РА 288°. Presumably, the accretion 
disk is also tilted toward this direction. Such a tilt is consis- 
tent with the EHT total intensity image shown in EHTC IV. 
Thus, if we measure angles counter-clockwise with respect 
to the presumed tilt direction in M87*, the polarized flux is 
strong between angles ~ +10° and —140°, with peak at —90* 
and —50° on April 5 and April 11. 

In our analytic model, the tilt and putative jet are toward 
the North. Thus, for a direct comparison of this model 
with the M87* image, we should rotate the calculated image 
clockwise by 72°. Alternatively, we could measure angles 
as offsets from the jet direction North. Thus, for a model to 
reproduce what is seen in M87*, it should have strong lin- 
early polarized flux between +10° from the jet, i.e., just to 
the left of North, and —140° from the jet, which is located in 
the lower-right quadrant. That is, the polarized flux should 
concentrate in the right half of the panels in plots such as 
Figs. 3-5 below, shading towards the upper right quadrant. 
As we will see, this is a fairly strong constraint. 

The second piece of information from the polarized im- 
age of M87* is that the polarization vectors show a twist- 
ing pattern that wraps around the black hole (EHTC VII; 
EHTC УШ). The twist is described quantitatively by the 
f» mode of the azimuthal decomposition of polarization de- 
scribed in Palumbo et al. (2020). The amplitude of £2 de- 
scribes the degree to which the EVPA obeys rotational sym- 
metry and scales linearly with fractional polarization, while 
the phase of 32 describes the twist angle between the EVPA 
and the local radial unit vector on the image. In the M87* 
image, the twist angle is fairly stable in the regions where the 
polarized flux is strong. With respect to the local radial direc- 
tion, the EVPA of the polarization vector is rotated clockwise 
by ~ 70°. This too is a strong constraint on models, as dis- 
cussed at length in EHTC VIII. 


3.1. Models with Pure Vertical Field 


Gravity Collaboration et al. (20182) reported observations 
of polarized flares in Sgr A* in near-IR, and showed that a 
model with a dominant vertical magnetic field can reproduce 
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Figure 3. Polarization patterns corresponding to models with a “vertical” magnetic field (non-zero В, in the fluid frame). In each case, the 
directions of the ticks indicate the orientation of the polarization É-vector around the ring as viewed on the sky. The lengths of the ticks are 
proportional to the polarized intensity. Top Left: Ring with a very large radius and no orbital velocity, so that neither velocity aberration nor 
lensing plays a role. Top Right: Large ring radius (i.e., no lensing), and fluid orbiting with a tangential velocity 8 = 0.3 in the clockwise 
direction (y = —90°). Bottom Left: Large ring radius (no lensing), and fluid flowing with velocity B = 0.3 radially inward (y = – 180°). Bottom 
Right: Ring with a small radius R = 6M, hence strong gravitational lensing, but with no fluid velocity, hence no aberration. 
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the observations. Motivated by this, we begin by studying the 
predictions of our toy model for a pure vertical field, oriented 
normal to the plane of the emitting ring. 

Figure 3 shows results from the analytical model for the 
case when В, = 1, B, = B, = О. It explores the two 
primary physical effects other than magnetic field direction 
that influence the observed polarization: (i) Doppler beam- 
ing and relativistic aberration caused by motion of the radiat- 
ing fluid, and (ii) gravitational lensing caused by the gravity 
of the black hole. The Top Left panel in Fig. 3 corresponds 
to a ring with a large radius (R = 10*) such that there is 
negligible gravitational lensing.We also set 8 = 0, thereby 
eliminating Doppler beaming and aberration. The only re- 
maining effect is the tilt of the ring, which causes the pure B; 
field in the ring frame to appear in projection on the sky as 
a vertically oriented (North-South) field. The polarized syn- 
chrotron emission from the ring has its EVPA perpendicular 
to the projected field, 1.е., in the East-West direction. The ob- 
served polarized intensity, which is indicated by the sizes of 
the polarization ticks in the plot, is uniform around the ring. 
In this figure and all others shown in this section, ticks are 
shown at 50 equally spaced positions in ф. 

The Top Right panel in Fig. 3 shows the effect of including 
an arbitrary relativistic velocity (8 = 0.3) for the fluid in the 
clockwise tangential direction (y = —90°), but still keeping 
a large radius, hence no gravitational deflection. In this case, 
there is a strong asymmetry in the polarized flux around the 
ring. However, the bright region of the ring is in the left half 
of the plot, exactly the opposite of what we require to explain 
MS87*. 'This contrary behavior is actually rather surprising. 
Given the direction of the tilt and the clockwise sense of ro- 
tation, the fluid in the right half of the plot has a component 
of its motion towards the observer, while the fluid on the left 
has a component away from the observer. Doppler beaming 
ought to favor the right side, yet we see the opposite. This 
paradoxical behavior is because of aberration, as we explain 
in sec. 4. 

The Bottom Left panel in Fig. 3 shows the effect of a pure 
inward radial velocity (y = —180°), again for a large ring 
radius. Once again, the bright region of the disk is on the 
wrong side compared to what is seen in M87*. It is also ex- 
actly the opposite of what we would expect from Doppler 
beaming, since the fluid in the upper half has a velocity com- 
ponent towards the observer, and ought to be bright. Once 
again, aberration is the explanation. 

Finally, the bottom right panel considers a ring at small 
radius (R = 6) such that gravitational deflection of light rays 
is important. For simplicity, we assume that there is no fluid 
velocity. In this case, the results are similar to the Bottom 
Left panel, and the strongest polarized flux is at the bottom, 
which does not match what is seen in M87*. 

We do not discuss the $2 phase of the polarization patterns 
for models with pure vertical field, except to note that in the 
regions where M87* has its strongest polarized flux (upper 
right), the sense of the EVPA twist seen in all the examples 
in Fig. 3 has the wrong sense. 


The conclusion from these examples is the following. If 
the polarized emission that we see in M87* at 230 GHz is 
from equatorial gas, and if the gas rotates in the clockwise 
direction, as EHTC V concluded, and/or flows radially in- 
ward, as is natural for accretion, then the magnetic field can- 
not be dominated by a pure vertical component. There must 
be substantial radial and tangential field components. 

Note that the observed ring in the Bottom Right panel in 
Fig. 3 has a radius slightly larger than the original ring radius 
R = 6. The ring is also shifted slightly upward relative to the 
origin. Both effects are the result of gravitational deflection, 
as we explain in sec. 4. The effect is seen only when R is 
small (gravity is strong), which is the case in this panel of 
Fig. 3, and in all the panels in Figs. 4, 5. 


3.2. Models with Pure Radial or Tangential Field 


We now turn our attention to models with magnetic field 
entirely in the equatorial plane, i.e., B; — 0, non-zero B, or 
By. We consider a ring with small radius (А = 6) and include 
relativistic fluid motion; thus, lensing, Doppler and aberra- 
tion are all included. Figure 4 shows four models, two with 
radial field (7 = 0°) and two with tangential field (7 = 90°). 
For each field configuration, we consider two velocity fields, 
either pure clockwise rotation (y = —90°) or pure radial infall 
(x = -180?). 

Three of the four panels in Fig. 4 have their strongest po- 
larized flux in the correct region of the ring (top and/or right) 
to match what is seen in M87*. Even the fourth (Top Right 
panel) has slightly stronger polarized flux at the top. The 
very different behavior of these models, compared to those 
in Fig. 3, is explained in detail in the next section. In brief, 
for models with magnetic field restricted to the equatorial 
plane, aberration induces the same sense of flux asymme- 
try as Doppler beaming and therefore enhances the effect 
of the latter, whereas in the pure В, models, aberration in- 
duces flux asymmetry with the opposite sign of that due to 
Doppler beaming, and in fact overwhelms the latter and re- 
verses the sign of what is observed. In this sense, equatorial 
field-dominated models are more promising for M87*. 

Considering the twist of the polarization pattern, as dis- 
cussed in EHTC VIII, a pure tangential field is ruled out be- 
cause the polarization ticks are predicted to be purely radial, 
which does not match M87. A pure radial field is also ruled 
out since it predicts polarization ticks entirely in the tangen- 
tial direction. However, these models come closer to what is 
seen in M87*. It would appear that models in which B, > Bg 
are most suitable. 


3.3. Models with Both Radial and Tangential Field 


Figure 5 shows four models in which both B, and B, are 
non-zero, and B; = 0. All the models have fluid with clock- 
wise rotation in the sky and radial infall, i.e., the angle y of 
the vector B is in the lower left quadrant. Since the radial and 
tangential magnetic field components in the inner regions of 
an accretion disk are likely oriented parallel to the motion of 
the fluid — the field is “combed out" by the flow — we sim- 
plify matters by assuming that the field is aligned with the 


POLARIZED IMAGE OF A RING ORBITING A BLACK HOLE 


0.= 20°, R = 6.0, B=0.30, y=-90°, B, Bg, Bz = 0.00,1.00,0.00 


05, 2 20°, R=6.0, В = 0.30, y=-180°, B,, Bg, Bz = 0.00,1.00,0.00 


8 8 
О] 
| cio EF : "Ur 
S "4 N А 
‘N P4 N a 
4 M p d 4 “чы Pd 
~ Шш TA — 
~ a — — 
24 _ ee ^ ___ Kem 
Sgt ы E "EN x 
^ m = z ph М 
-2 g i -2 p ae 
S x e D 
е ^ N E Pd N 
/ N / N 
^4, De // GS 
i Bb r E 39 
-8 -8 
-8 -6 E E 0 2 4 6 -8 -6 -4 =2 0 2 4 6 
x(M) x(M) 
8, = 20*, R - 6.0, 8 = 0.30, у =-90°, В, By, B; = 1.00,0.00,0.00 8, = 20*, R - 6.0, 8 = 0.30, = -180°, Br, By, B; = 1.00,0.00,0.00 
8 8 
P ER 
6 d 6 
7 
"A 
4 vA 4 
/ 
2 / 2 
A l £ 
È XE: | 
| 
\ \ / 
5 \ = \ / 
\ \ / 
E \ P \ / 
SE NS 7 
-6 эмас 2 -6 Loc ш 
-8 -8 
-8 Z6 En E 0 2 4 6 -8 Z6 4 29 0 2 4 6 
x(M) x(M) 


Figure 4. Polarization patterns for models with magnetic field in the equatorial plane. Top Left: Azimuthal field ( = 90°) with azimuthal 
clockwise velocity (y = —90°). Top Right: Azimuthal field (7 = 90°) with radial inward velocity (y = —180°). Bottom Left: Radial field 
(7 = 0°) with azimuthal clockwise velocity (y = —90°). Bottom Right: Radial field (7 = 0°) with radial inward velocity (y = – 180°). 
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velocity. Specifically, we choose 


Pure Beg: =X Or N=VtT. (31) 


For the specific case of a purely equatorial field, we can 
choose either of the two values of 7 indicated above. The 
two choices correspond to oppositely oriented directions of 
the magnetic field lines; this ambiguity has no effect on the 
linear polarized emission. As we discuss in subsection 3.5, 
we need to be more careful about the choice of 7 when we 
have both vertical and equatorial field components. 

In Figure 5, the model in the Top Left panel has tangen- 
tial velocity larger than radial velocity, and correspondingly 
By > B,. In the Top Right panel, the radial and tangen- 
tial components are equal, while in the lower two panels the 
radial components of velocity and magnetic field are larger 
than the respective tangential components. All four models 
have flux asymmetry that qualitatively matches M87*. All 
four models also have polarization patterns with the same 
sense of twist, or sign of 6) phase, as observed in M87*. 
Among the four models, the ones in the bottom row come 
closest to M87*. 


3.4. Models with R = 4.5 M and Varying Inclination 


We round out the discussion of examples by considering 
models with a smaller emission radius, R = 4.5, which is 
better matched to M87*, and exploring the effect of varying 
the tilt angle 6. Figure 6 shows models with y = —150°, 
7] = x + л = 30°, and four choices of 04: 20°, 40°, 60°, and 
80°. 

The top left panel has 0, = 20° and is designed to resem- 
ble M87*. The polarized intensity asymmetry (relative to the 
direction of the jet), as well as the twist of the EVPA pattern, 
are similar to the EHT observations described in EHTC VII 
and EHTC VIII. This same model is shown again in Fig- 
ure 9 with the polarization pattern rotated counter-clockwise 
by 288° to match the jet orentation in M87*, and with the 
emitting fluid spread out in radius with an exponential profile 
with scale width 2M (see subsection 5.1 for details), instead 
of the infinitely thin emitting ring assumed here. 

The remaining panels in Figure 6 show the effect of in- 
creasing the tilt angle 0. The Doppler asymmetry in the po- 
larized intensity increases rapidly since the fluid motion has 
a larger component parallel to the line-of-sight. The orienta- 
tion of the asymmetry (bright on the right, dim on the left) as 
well as the twist of the polarization pattern qualitatively re- 
semble what is seen in the 0, = 20° model. The ring appears 
increasingly flattened as 0, increases, but it also acquires an 
additional asymmetry such that, by 0, = 80° it looks more 
like a semi-circle than an ellipse. This is because of extreme 
lensing of radiation emitted from the far side of the ring. As 
in the previous Figures, ticks are equally spaced in ¢; the 
large gaps on the north side of the Ө = 80° image indicate the 
relative stretching between ф and ф at high inclination. 


3.5. Models with All Field Components 


We finally discuss models in which all three components 
of the magnetic field are non-zero. In this general case, 


we need to be careful about the geometry of the magnetic 
field. In a three-dimensional accretion flow in which mag- 
netic field lines penetrate the disk from one side to the other, 
as for instance in a magnetically arrested disk (MAD) field 
geometry (Narayan et al. 2003; Igumenshchev et al. 2003; 
Tchekhovskoy et al. 2011; Bisnovatyi-Kogan 2019), one ex- 
pects a reflection antisymmetry іп B,, about the midplane. 
That 15, B, and В, would flip sign when crossing the mid- 
plane, whereas В, would retain the same sign on the two 
sides. Let us assume, without loss of generality, that B. is 
positive, i.e., the z-component of the magnetic field line is 
pointed towards the observer, and let us also take B, to be 
positive. If the magnetic field is dragged and aligned with 
the flow, as we assumed in the previous two subsections, the 
field angle ту and the flow velocity angle y must be related as 
follows on the two sides of the disk, 


z > 0 (near side): n=y+7, 
z <0 (far side): n=y, (32) 


where “near side” means the side of the disk facing the ob- 
server. 

In the absence of Faraday rotation effects, the above anti- 
symmetry affects emission only by changing the relative sign 
between Beq and B,, hence it is not relevant if either Beg or 
В, is zero. However, when both Beq and В, are non-zero, опе 
should separately compute the polarized image produced by 
the near side and far side of the disk and add the resulting 
Stokes parameters. 

If Faraday effects internal to the flow are strong enough to 
depolarize the emission from the far side, the polarized image 
seen by the observer will be dominated by the near side. The 
simulations considered in EHTC УШ, for instance, generally 
show large internal Faraday depths. In such cases, we need 
compute only a single image from the near side of the disk, 
setting n=yt+7. 

We do not show examples of models with both vertical and 
equatorial field since the parameter space is large. 


3.6. Numerical Geodesics and Effect of Spin 


A general Beloborodov-like analytic approximation for the 
emission angle of photons from equatorial matter around a 
spinning black hole is not known. However, it is possible to 
solve analytically for the observed polarization once the pho- 
ton’s arrival coordinates on the image are determined from 
a numerical solution to the geodesic equation; this relation 
can be explicitly expressed in terms of real elliptic integrals 
(Gralla & Lupsasca 2020a,b, see also Li et al. 2005; Gates 
et al. 2020 for a calculation of images of an orbiting emit- 
ter in this formalism). For a spinning black hole, we gener- 
alize the P-frame to the “zero-angular-momentum-observer” 
(ZAMO) frame, and then consider a boost B as in (8) into 
the corresponding F-frame. The semi-analytic result for the 
polarized image of such a boosted fluid orbiting a spinning 
black hole is presented in Figure 7, in which changing spin 
is plotted by color. The inner and outer ring in the first two 
panels correspond to emission radii of R = 4.5 and R = 6, 
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Figure 5. Polarization patterns for four models that include both radial and azimuthal components of velocity and magnetic field. The models 
correspond to y = —120° (Top Left), y = —135° (Top Right), y = – 150° (Bottom Left), y = —165? (Bottom Right), each with magnetic field 
trailing opposite to the velocity (7 = y + 180°). The two models in the bottom row come closest to reproducing the polarization pattern seen in 
M87*. 
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respectively. The results of the Beloborodov approximation 
are overlaid with black dashed lines and coincide with the 
low spin semianalytic solution from Kerr. The first and sec- 
ond panels of Fig. 7 generalize the scenarios from the bottom 
right panel of Fig. 3 and the upper left panel of Fig. 5, respec- 
tively. The small panels zoom in on one set of ticks from the 
second panel. 

Figure 7 illustrates that for the idealized case of purely ge- 
ometric and relativistic effects that we consider here, black 
hole spin has only a small effect on the observed EVPA 
and can be reasonably neglected for the purposes of the toy 
model. It also shows that the Beloborov approximation is 
fairly accurate even at radii as small as R = 4.5. The effects 
of spin on observed polarization become more pronounced 
at very small radius and high observer inclination, neither of 
which are considered in this paper but will be the subject of 
future work. 


3.7. Generalizations 


Although the examples presented in this paper are re- 
stricted to axisymmetric models with emission limited to a 
single radius, the underlying model is more general. The pri- 
mary result of the analysis presented in sec. 2 is an analytical 
method to map emission properties at a given (А, ф) in the 
emitting ring to the properties of the observed radiation in 
the sky plane. This transformation can be easily applied to 
models with non-axisymmetric emission, as well as to radi- 
ally extended sources. In such models, IBI would be a func- 
tion of location and this would need to be included in the 
calculations. Other quantities like the electron temperature 
and number density that affect the emissivity could also vary 
with position and will need to be accounted for. 

Two other approximations in the model, both made in the 
interests of simplicity, deserve discussion: (1) We restricted 
the emitting gas to lie in a single equatorial plane. (2) We 
took the velocity to lie entirely within the same plane (though 
we did allow for a general magnetic field). Both limitations 
can be eliminated. 

The Beloborodov approximation can be applied at any 
emission location (R, ф, 2), not just at equatorial locations. 
For non-equatorial locations, the geometry of the Geodesic 
Frame and the computation of a (Fig. 1) will differ. This will 
modify the result for the components of Eis If a given null 
geodesic has contributions from several emission regions at 
different heights z from the equatorial plane, one could com- 
pute their individual contributions to the Stokes parameters 
and add the contributions incoherently. 

Similarly, an off-plane velocity component will modify the 
Lorentz transformation coefficients between the P-Frame and 
the F-Frame, and will alter the geometrical factor that enters 
the path length calculation. The distinction between “ver- 
tical” and “in-plane” magnetic field components would be- 
come less clear, but this is merely a matter of definition. 

The model discussed in this paper has been derived for 
a non-spinning (Schwarzschild) black hole. However, as 
shown in subsection 3.6, and as discussed also in Gravity 
Collaboration et al. (2020) and EHTC VIII, black hole spin 


has very little effect on the polarized image, at least for the 
low inclination angles considered so far. 

Finally, the analysis here is focused on optically thin syn- 
chrotron emission for which the polarization four-vector f" 
is given by equation (16) and the electric field is normalized 
as in equation (25). For optically thick emission from a thin 
accretion disk, other prescriptions will need to be substituted, 
e.g., Li et al. (2009) discuss polarization of X-rays emitted by 
the scattering atmosphere above a black hole X-ray binary 
disk. Except for this change, the rest of the analysis should 
remain the same. 


4. ANALYTICAL UNDERSTANDING OF THE 
RESULTS 


By Taylor-expanding the expressions given in sec. 2 in 
suitably chosen “small” quantities, and keeping terms up to 
second order, we can obtain useful analytical approximations 
for various observables. This provides a physical understand- 
ing of the results shown in sec. 3. 

In the present context of trying to understand M87* and 
Sgr A*, we have three small quantities, 2/R « 1/3 (lensing), 
В = 1/3 (Doppler and aberration), sin@, ~ 1/3 (ring tilt), 
where the numerical values correspond to the models shown 
in sec. 3. We treat all three quantities on an equal footing in 
the series expansions we carry out. The full results, with all 
terms up to quadratic order, are listed in Appendix D. The 
reason for going up to quadratic order is explained below. 
Here we use the series expansion of the equations to interpret 
the numerical results presented in sec. 3. 


4.1. Shape of the Observed Ring 


We begin with the shape of the ring as observed on the sky. 
To quadratic order, the result is 


x =(R+1)cosy (33) 
+ = COS Q + sin 6, sin2y — 2 sin? 0, sin? ф сов |, 
2R 2 
y=(R+1)sing (34) 
+ ВЕРИ Ga’ dno Р 
2R 2 


The first term in each expression gives the answer up to lin- 
ear order, and the remaining terms inside the square brack- 
ets correspond to quadratic order. Up to linear order we see 
that the observed ring is circular, but with an apparent radius 
larger by unity (i.e., GM/c?) than the radius of the source 
ring. The radial “expansion” of the observed ring is caused 
by gravitational deflection (lensing) of geodesics. As shown 
in Fig. 1, lensing causes the geodesic to curve around the 
black hole such that the impact parameter is larger than the 
naive straight-line estimate А sin y. 


5 In the case of M87*, observations of the radio jet suggest a tilt б, ~ 17° 
(Walker et al. 2018), and in the case of Sgr A*, Gravity Collaboration et al. 
(2018a) estimate 6, < 30° based on the polarization signatures of infrared 
flares. 
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Figure 6. Polarization patterns for four models with equatorial magnetic field and emission radius R = 4.5, viewed at different inclination 
angles. Top left: 6, = 20°. Top right: 0, = 40°. Bottom left: 0, = 60°. Bottom right: 0, = 80°. All the models have velocity angle y = – 150°, 
and magnetic field trailing opposite to the velocity (7 = y + 180°). The model in the top left, rotated counter-clockwise by 288° and with 
emission spread over a finite range of radii, is shown in Figure 9 as a toy model of M87*. 
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Figure 7. The effects of spin on the observed polarization pattern. Each of the two main panels displays a different configuration of magnetized 
fluid. The first panel corresponds to the bottom right panel of Fig. 3 and the second panel corresponds to the top left panel of Fig. 5. Both panels 
show an inclination of 20° and negative spin (i.e., clockwise rotation on the image). The inner and outer rings of polarization ticks correspond to 


emission from R = 4.5 and R = 6, respectively. The color bar shows increasing spin from a = 0 to |a| = 1, and the Beloborodov approximation 
for Schwarzschild is shown in black overlaid dashes. The two small panels display a zoom-in of one set of ticks at R = 4.5 (lower) and R = 6 


(upper). 


Among the quadratic terms in equations (33) and (34), the 
terms proportional to 1 /R are second-order corrections to the 
ring radius, and the sin? 6, terms describe the flattening of the 
observed ring because of tilt. The latter is simple geometry: 
a tilted circular ring appears elliptical in shape, with a minor 
axis radius equal to cos 6, = 1 — (1/2) sin? 6, times the orig- 
inal ring radius. The sin 6, terms describe the effect of tilt 
on lensing. Geodesics reaching the observer from the upper 
half of the ring (0 < ф < л) travel a longer distance near the 
black hole and suffer more deflection (this is the case shown 
schemaically in Fig. 1), while geodesics from the lower half 
(x < ф < 2л) experience less deflection. This causes an up- 
ward shift of the observed ring, i.e., a net positive bias in y. 
The shift is of the order of sin 6, in units of GM/c?. The shift 
is seen in all the models in sec. 3 that have a smallish radius 
(R = 6, Lower Right panel in Fig. 3, and all panels in Figs. 4, 
9,0). 


4.2. Doppler Factor and sin С 


Expanding up to second order, we find for the Doppler fac- 
tor ó, 


1 
б= ( - 5] 
R 
2 1 2, 

-|£ + 27720 E cosy + Bind sin(y + 9) ‚ (35) 
where the second order terms are shown on the second line 
inside square brackets. The linear order term —1/R describes 
deboosting of the observed intensity by gravitational red- 
shift, and the first three second-order terms describe various 
other deboosting effects such as second-order Doppler. Since 


cos y is negative for radial infall, all three terms have a posi- 
tive magnitude for the inflowing models we have considered, 
causing uniform dimming all around the ring. 

Azimuthal modulation of the intensity from relativistic 
beaming is described by the final term, £ sin 6, sin(xy + Ф), 
and this is the only term that varies as a function of o. The 
fact that this important effect appears only at second order is 
a major reason for expanding the equations up to quadratic 
order rather than stopping at linear. Why is it second order? 
It is because azimuthal modulation from Doppler beaming 
requires both tilt and fluid velocity, each of which is treated 
as a small quantity in our analysis. 

Doppler beaming causes an increase in the observed polar- 
ized intensity when sin(y + 4) is negative, with the maximum 
boost occurring when y + y = —90°. For pure clockwise 
rotation (y = —90°), the maximum boost is at y = 0. This 
is natural since, for a ring tilted towards the North, the fluid 
at o = 0 has the largest velocity component towards the ob- 
server and hence produces the most Doppler-boosted radia- 
tion. For pure radial infall (ү = —180°), the maximum boost 
is at o = 90°, again because the fluid there has the maximum 
velocity towards the observer. Since we consider models that 
lie between these two extremes, we expect the polarized in- 


For the models considered in sec. 3, where each of the three small quan- 
tities is ~ 1/3, one expects second-order terms to be of order 10% of the 
leading-order terms. However, many second-order terms come with large 
coefficients, e.g., intensity is proportional to ó^ so Doppler boost goes like 
—4ß sin, sin(x + y). Hence the second-order contributions are often not 
small. The analysis in this section should thus be used only for qualitative 
understanding. For accurate results, it is necessary to evaluate numerically 
the full equations given in sec. 2. 


POLARIZED IMAGE OF A RING ORBITING A BLAcK Нов 17 


tensity to be maximum somewhere in the top right quadrant, 
0 < ф < 90° (for a tilt to the North). This agrees with what 
is observed in M87* (once we allow for the different tilt/jet 
direction). Surprisingly, it is not true for the models shown 
in Fig. 3. To understand the reason for this discrepancy, we 
need to consider a second effect. 

From equation (15), the observed polarized intensity de- 
pends on the Doppler factor 6 as well as the path length /, 
and the angle Z between the photon wave-vector Ер in the 
fluid frame and the local magnetic field B. For small tilt an- 
gles, the variation in the path length is small and not very 
important. We ignore it in the discussion below. The angle 4, 
however, is crucial since synchrotron emission is maximum 
when Ер апа В аге orthogonal to each other (4 = +л/2) 
and vanishes when they are parallel (¢ = 0, л). Appendix Р 
evaluates 182 sin? £ up to quadratic order. We consider іп the 
following subsections the effect of various terms in the series 
expansion. 


4.3. Models with Pure Vertical Field 


We begin by considering a model with pure B, and con- 
sider the non-zero terms in I8? sin? С: 


B, Finite, Beg = 0: 
4p 


4 4 
|BI sin? 4 = ^R sin 6, sin ф + m + sin? 0, — R COS Y 
+ 26 sin 6, sin(y + ф) + B^--- B2. (36) 


There are several interesting effects here. First, we have only 
second-order terms, no zeroth- or first-order terms (this is 
another reason for going up to second order in the analy- 
sis). It suggests that the observed flux should be strongly 
suppressed. This is not surprising since the emission towards 
the observer goes as sin? Z ~ sin? 6,, which is small for mod- 
els with small tilt. The lack of zeroth- and first-order terms 
also means that the importance of the second-order quantities 
in equation (36) is enhanced. 

Consider first the term —(4/R) sin Ө, sin y, which describes 
the combined effect of lensing (4/R) and tilt (sin 9,). Fig- 
ure 1 shows the origin of this term. In the absence of lensing, 
a geodesic travels on a straight line to the observer and hence 
subtends an angle 6, to the (vertical) magnetic field. When 
gravitational ray deflection is included, the angle at the emis- 
sion point is modified. For a point on the North or upper half 
of the ring (the case shown in Fig. 1), the deflection is such 
that the photon wave-vector becomes more nearly parallel to 
the z-axis, i.e., more parallel to the magnetic field. Thus Z is 
reduced, and this causes the emissivity to go down. The de- 
crease is largest when y = 90°, as indeed we find in equation 
(36). If we consider instead a point on the South or lower 
half of the ring, e.g., o = —90°, the gravitational deflection 
works in the opposite sense and causes Z to increase, and 
the emissivity to correspondingly increase. The net result is 
an asymmetry in the polarized flux around the ring such that 
the maximum flux is in the South and the minimum is in the 
North, precisely as seen in the Bottom Right panel in Fig. 3. 


Consider next the term 26 sin 6, sin(y + ф), which corre- 
sponds to the combined effect of tilt and relativistic motion. 
Here the relevant effect is aberration. Because of the mo- 
tion of the fluid, the orientation of the wave-vector Ёр) in the 


D . . . . . p . 
fluid frame is different from its orientation kp) in the P-frame. 
The aberration effect is such that fluid that is moving towards 


the observer has Ер rotated closer to е z-axis in the fluid 
frame, i.e., more nearly parallel to B, while fluid that is mov- 
ing away from the observer has the tilt of Ер with respect 


to B increased. The former fluid element thus emits less and 
the latter more in the direction of the observer. This cancels 
the effect of Doppler beaming. Actually, since the constant 
q-independent terms in equation (36) are of the same order as 
the modulation term sin(y + ¢) (note that 28 sin 0, is almost 
equal to 4/R? + sin? 6, +), the cancellation tends to be quite 
pronounced when y + ф ~ —90°. The net effect is that aber- 
ration overwhelms Doppler beaming and gives the patterns 
seen in the Top Right and Bottom Left panels in Fig. 3. 


4.4. Models with Pure Equatorial Field 


When we consider models with pure equatorial field (Beg 
finite, B, = 0), the situation is quite different. Focusing on 
IB? sin? £, we find 


Beg Finite, B, = 0, п = y t7: 


IBP sin? z & B2, + [-28sin&, sin(y + )--+] В, (37) 
where we have written only one of the second-order terms. 
As in sec. 3, we have simplified matters by assuming that 
the magnetic field is oriented anti-parallel with the velocity: 
n-2xy-m. 

The first thing to note is that in the case of an equatorial 
field there is a non-vanishing zero-order term. For small tilt, a 
magnetic field in the equatorial plane is almost orthogonal to 
the photon wave-vector, hence synchrotron emissivity in the 
direction of the observer is nearly maximum. Correspond- 
ingly, the second-order terms are less important. Moreover, 
the second order term in equation 37 appears with the same 
sign as the corresponding term in ó (eq. 35), and the opposite 
sign as in equation (36). The reason is simple. When aberra- 
tion tilts the wavevector closer to the z-axis, the wavevector 
becomes more nearly orthogonal to B, and hence the emis- 
sivity increases. Thus in equatorial field models, the second- 
order terms in 182 sin? 4 cooperate with and enhance the ef- 
fect of Doppler beaming, as seen in the panels in Figs 4 and 
5. As an aside, when both B, and В, are non-zero, and if we 
assume as before that 7 = y +z, then there is a first order term 
—2sin 0, sin(1] + p) ВВ, which again has the same sign as 
the corresponding term in б. 


4.5. Twist of the Polarization Pattern 


We now briefly discuss the twist of the polarization pattern 
around the ring. When the field is purely in the equatorial 
plane, the results are transparent. To zeroth order, the electric 
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field in the sky plane is given by 


Exons = = Sing B, — cos Ву = = sinf + Ф) Beg, 
Ез = cosp B, — sin g By = cost) + Ф) Beg (38) 


That is, the electric field is oriented perpendicular to the pro- 
jected magnetic field, as one would expect. 

Instead of considering the electric field, one could consider 
the Stokes parameters Q and U and look at their Fourier 
coefficients Bm (Palumbo et al. 2020), as described in Ap- 
pendix D. The most useful coefficient is 62, whose complex 
phase directly gives the orientation of the twist. If the electric 
field is radial, the phase of 32 is zero, if it is rotated clockwise 
from radial by 45°, the phase is —90°, and if the electric field 
is tangential, the phase is –180°. The ЕНТ observations of 
M87* give a phase ~ —130? = +230°. From Appendix D, 
the leading order term in [82 in the case of a pure equatorial 
magnetic field is 

pase Be (39) 


The phase of this quantity will match the phase observed in 
М87* if 7 ~ 25°. Hence, the magnetic field must be mostly 
radial. 

When B,, = 0 and we have a purely vertical field, the phase 
of 82 is determined by the coefficient of B), which consists 
entirely of second-order terms: 


Ве =0: В = lox + AB i = een) B? . (40) 


к К 


If lensing is unimportant, i.e., R is large, then 82 dominates 
and the phase of 82 is determined by the orientation angle y 
of the fluid velocity. For a radial velocity (y = л), the phase 
of f is л, i.e., the polarization vectors should be tangentially 
oriented. This is indeed seen in the brightest part of the ring 
in the Bottom Left panel in Fig. 3. Similarly, for a tangential 
velocity (y = —7/2), the phase of B; = 0 and the polarization 
ticks should be radial, as seen in the Top Right panel of Fig. 3. 
Finally, if there is no velocity but we consider strong lensing 
(small А), then equation (40) shows that 62 has phase = л and 
the polarization should be tangential, as in the Bottom Right 
panel. 


5. COMPARISON TO OBSERVATIONS 


Our ring model provides a convenient framework for di- 
rect comparison with a variety of polarimetric observations 
of near-horizon emission. We now discuss two specific cases 
of particular interest: polarimetric imaging with the EHT and 
infrared flares of Sgr A*. 


5.1. Comparison to the M87 Polarized Image 


Recent EHT observations produced polarized images of 
M87* (EHTC VID. As reported in the one-zone model com- 
parisons performed in EHTC V and EHTC УШ, the bright- 
ness, angular size, and expectation of significant Faraday ef- 
fects coarsely constrain the magnetic field strength B, elec- 
tron number density ne, and electron temperature Т, in the 


flow imaged by the EHT. The EHTC VIII results suggest that 
B € 30G, 10^ <n, < 10cm, and 10? < T, < 1.2x10!! K. 
The reconstructed images in EHTC VII were compared to 
general relativistic magnetohydrodynamic (GRMHD) sim- 
ulations to identify a space of favored model parameters 
(EHTC УШ). We will now explore whether our ring model 
can reproduce the polarization structure in these favored 
GRMHD simulations and in EHT images of M87*. 

For the GRMHD comparison, we first perform an az- 
imuthal and temporal averaging in the fluid domain to ap- 
proximate a stationary axisymmetric flow. In the fluid frame, 
the magnetic field in each cell is decomposed in Cartesian 
Kerr-Schild coordinates, which are then recast into cylindri- 
cal coordinates and then azimuthally averaged. These az- 
imuthally averaged magnetic field decompositions are then 
further averaged over time between 7500 < t/(GM/c?) < 
10000 (the final quarter of these simulations). We then sam- 
ple values of the fluid velocity and magnetic field vectors 
from the averaged simulations and use these values to gen- 
erate ring models at Ө) = 17°. To avoid sampling near where 
the tangential and radial field directions tend to abruptly flip 
sign, we use z — 1M, just above the midplane. We use 
R = 4.5M, corresponding to the apparent lensed size of the 
emission ring in EHT images of M87* (see the later discus- 
sion of the observed image). To create an image from the 
one-dimensional ring model, we adopt a radial profile that 
decays symmetrically in R about R = 4.5 as an exponen- 
tial with a scale width of 2M (EHT images only constrain 
this width to be «5M; EHTC VI). We take a pixel-wise 
fractional polarization |m| of 0.7 before blurring in the ring 
model. Finally, we convolve both the ring model image and 
the GRMHD image with a 20 uas Gaussian kernel. 

Using this approach, Figure 8 compares four favored 
GRMHD models to the corresponding ring models. In each 
case, the ring model reproduces the sense of EVPA twist 
and relative polarized intensity of the averaged and blurred 
GRMHD image, although discrepancies in arg(82) suggest 
contributions from emission away from the midplane or from 
other effects that are not included in the ring model (e.g., 
black hole spin or Faraday effects). Тһе Riow and Rhigh pa- 
rameters adapted from Moscibrodzka et al. (2016) for use in 
EHTC V tune the ratio of electron to ion temperatures de- 
pending on the magnetic energy density of the plasma; large 
values of Rhign tend to produce significant emission far from 
the midplane, particularly in SANE models. Also, Faraday 
effects in MAD models can produce significant coherent 
rotation of the EVPA and, hence, in arg(82) (EHTC УШ). 

Figure 9 compares a representative ring model to the “соп- 
sensus" ЕНТ polarimetric image for 11 April, 2017 (1.е., the 
method-averaged image, see EHTC VII). The ring model pa- 
rameters are chosen based on the observed image and a priori 
expectations for M87*. For simplicity, we take B. = 0, al- 
though non-zero values of B,/Beq over a modest range also 
give similar results. We use y = —150°, to roughly match 
the observed 62 for M87* (see subsection 4.5). We take 
К = d/(20,) — 1 = 4.5 (subsection 4.1 explains the —1 fac- 
tor) where d = 42 uas is the observed ring diameter and 
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Figure 8. Comparison of GRMHD simulations to images of the ring model for simulation parameters favored in EHTC VIII. The left three 
columns show random snapshots, time averaged images, and blurred time averages of each GRMHD simulation; the right column shows the 
image generated by the simple ring model when evaluated for magnetic field and fluid velocity values taken from the simulations at R = 4.5 after 
azimuthal and temporal averaging. Ticks show polarization magnitude and position angle where total intensity exceeds 5% of the maximum. 
Grayscale shows total intensity in linear scale (directly proportional to polarization magnitude for the ring model). The total intensity and 
polarization magnitude are separately normalized in each panel. Panels show the average fractional polarization weighted by total intensity 
at bottom left; note that the GRMHD images are heavily depolarized, whereas the ring model images are not. The ring model and averaged 
images show the argument of the 8. PWP mode at top left. 
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Figure 9. Comparison of the EHT polarimetric image of M87* on 2017 April 11 (left) with a representative ring model (right). Ticks show 
polarization fraction (color), magnitude (length), and position angle (direction); grayscale is identical for the two panels and shows total 
intensity of the EHT image of M87*. Ticks are only plotted where the M87* polarization exceeds 296 of the maximum intensity. All images 
are shown after convolution with a circular beam of FWHM 23 uas (shown in the left panel). As in Figure 8, the total intensity and polarization 
are individually normalized for each panel. The ring model has clockwise rotation with radial inflow, corresponding to the top left model in 
Figure 6 after counterclockwise rotation by 288°. For complete model details, see subsection 5.1. The fractional polarization of the resolved 
ring model is set to 7096; the fractional polarization is reduced only through beam depolarization. Even after blurring, the ring model has 
significantly higher fractional polarization than the M87* image, although the relative variation in fractional polarization is similar across both 
images. 
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0, = 3.8 pas is the angular gravitational radius (EHTC VI). 
We use B = 0.4, which is comparable to the equatorial veloc- 
ity seen in GRMHD simulations (see Ricarte et al. 2020). We 
use ĝo = 20° to match the jet inclination of M87*. Thus, this 
model has a modestly relativistic fluid with clockwise rota- 
tion and predominantly radial infall. This model corresponds 
to the top left panel of Figure 6 after rotation to match the jet 
position angle of M87*, 288°. As with the GRMHD compar- 
ison, the ring model is evaluated over an exponential profile 
with a scale width of 2 M centered at R = 4.5 M. The result- 
ing ring model image is broadly consistent with the polariza- 
tion morphology of the EHT image. 

Although the qualitative agreement in Figure 9 is encour- 
aging, our simple ring model fundamentally fails to repro- 
duce all the features in the M87* image. Namely, our sim- 
plest model would produce a high fractional polarization 
(260%), while the M87* image has a low resolved fractional 
polarization € 2096. This suggests that significant depolariza- 
tion from internal Faraday effects are essential when model- 
ing and interpreting the M87* image. Nevertheless, the suc- 
cess of the ring model in reproducing the structure of some 
GRMHD images that have significant Faraday effects is en- 
couraging for the prospects of physical inference from this 
simple model. 

One possibility for using our model for a more complex 
emission scenario is to combine multiple ring models that 
correspond to different emission regions. Specifically, the 
assumption 7 = y + л corresponds to emission sourced by 
entrained magnetic field lines on the near side of the accre- 
tion flow (see subsection 3.5). The far side of the flow would 
instead have n = y, flipping Boq: Ignoring that contribution 
is equivalent to assuming that Faraday depolarization effects 
in the midplane are strong, so that the far-side emission is 
fully depolarized (as indicated in many models considered in 
EHTC VIII; see Ricarte et al. 2020). Our ring model could 
also be adapted to the case of weak Faraday rotation in the 
midplane; the resulting image would be the sum of two ring 
models, one with 7 = y and the other with 7 = y + л. Both 
cases would reduce the image polarization substantially and 
may give better agreement with the M87* image, but we de- 
fer a full analysis to a future paper. 


5.2. Comparison to Sgr A* Polarization 


The polarization of Sgr A* shows continuous variability in 
the submillimeter (Marrone et al. 2006; Johnson et al. 2015; 
Bower et al. 2018) and also shows rapid variability during 
near-infrared (NIR) flares (Eckart et al. 2006; Trippe et al. 
2007; Zamaninasab et al. 2010; Gravity Collaboration et al. 
2018b). The variability often appears as "loops" in Stokes Q- 
U, and is frequently attributed to localized emission from an 
orbiting “hotspot” (Broderick & Loeb 2005, 2006; Fish et al. 
2009). For the case of NIR flares, Faraday effects, absorp- 
tion, and background emission are insignificant, so we can 
directly compare observed values of polarization and cen- 
troid motion with a simulated hotspot-only model. 

Figure 10 shows a representative example. In this figure, 
we compute the hotspot polarized flux in the (Q, U) plane 
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Figure 10. Polarization signatures for a vertically magnetized 
hotspot on a circular, relativistic Keplerian orbit. Each curve shows 
the polarized flux for a full orbit. Different curves correspond to 
varying the hotspot radius (top) and viewing inclination (bottom). 
Note that we use radio astronomy conventions for Q and U here, 
distinct from those in Equation D.10 by an overall sign. 


over a full period for a set of orbits with varying emission ra- 
dius and inclination. We hold the underlying magnetic field 
structure to be vertical and constant, and adopt a relativistic 
Keplerian velocity for the hotspot: 8 = 1/ Vr — 2. Our results 
are similar to previous studies with fully numerical calcula- 
tions (see, e.g. Fish et al. 2009; Gravity Collaboration et al. 
20182, 2020); lensing and aberration compress the image of 
azimuthal evolution of polarization on one side of the flow 
and expand it on the other. In the formalism of azimuthal 
Fourier modes on the ring (Palumbo et al. 2020), power is 
shifted from the m = 2 mode to ће m = 1 mode. 
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We have developed an analytical method for computing the 
polarized image of a synchrotron-emitting fluid ring orbiting 
a Schwarzchild black hole. Given simple assumptions for the 
magnetic field geometry and fluid velocity, this model allows 
us to generate predictions of EVPA and relative polarized 
intensity as a polar function in the observed image at arbi- 
trary viewing inclination. We explored the main features of 
the model through a number of representative examples and 
by further expansion in the inverse emission radius (lensing), 
fluid velocity (Doppler and aberration), and observer inclina- 
tion (ring tilt). These reveal how the various physical effects 
influence the polarized image. 

In its simplest form, the fractional polarization of our 
model is significantly higher than that seen in EHT images of 
M87* (EHTC VII). This may indicate significant sub-beam 
depolarization, potentially from strong internal Faraday ef- 
fects (EHTC УШ). If so, observations at higher frequencies, 
where Faraday effects are suppressed, may show significantly 
higher image polarizations, while observations at lower fre- 
quencies are expected to show a heavily depolarized “core.” 

Our polarized ring model provides intuition and insights 
about how a black hole’s accretion flow and spacetime com- 
bine to produce a polarized image. It also provides a pathway 
to constrain these physical properties through direct compar- 
isons with data and images from the EHT, GRAVITY, and 
future X-ray polarimetry studies. Extensions such as non- 
axisymmetric structure and non-equatorial emission will pro- 
vide an expanded class of geometrical models to complement 
the growing library of GRMHD simulations (EHTC V). The 
inclusion of black hole spin will be necessary for rigorous un- 
derstanding of M87* polarization, particularly if emission at 
small radii is significant. Further studies which examine the 
capability of the model in matching snapshots of GRMHD 
simulations with similar magnetic field and flow conditions 
will elucidate how readily field geometries may be directly 
inferred from polarized images. 
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Figure 11. Testing the accuracy of the Beloborodov approximation. The left panel shows fractional error in p = yx? + y? as a function of 
y for р = 3, 5, 7, and 9. Yellow ranges denote values of y relevant for observer inclinations 0, = 20° and 60°. The center and right panels 
show the image coordinates for rings with emission radius R = 2 (red), 4 (green), 6 (blue), and 8 (cyan) viewed at inclinations of 20° and 
60°, respectively. For each ring, the solid line shows the exact calculation, while the dotted line shows the Beloborodov approximation (see 


Equation C.4). The black dotted line shows the critical curve, p = р, = V27. 
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APPENDIX 
A. ACCURACY OF THE BELOBORODOV APPROXIMATION 


The model developed in section 2 relies on the approximate formula Equation 4 derived by Beloborodov (2002). This approx- 
imation provides an estimate for œ (and, equivalently, for р; Equation C.4) for given emission coordinates К and ф. We now 
quantify the accuracy of this approximation. 

Emission from the equatorial plane arriving at a given observer inclination angle 0 € 0, < 7/2 will sweep through y € 
[1/2 + 0,] as the azimuthal angle ø varies (see Equation 1). In particular, all emission from a face-on disk has y = 7/2, while 
emission from an edge-on disk samples angles 0 < y < л. As the left panel in Figure 11 shows, the error in the Beloborodov 
approximation increases with y. In the context of the ring model, the approximation is most accurate at small inclinations. For 
0, = 17° for example (relevant for M87*), the approximation for p has a fractional error smaller than 2% for all values of R. This 
error decreases rapidly as p grows; e.g., for p = 9, the fractional error in p is smaller than 0.03%. In general, for emission on the 
side of the accretion disk closer to the observer (i.e., л < ф < 2л, у < 7/2), the approximation for p will have fractional error 
smaller than 0.6% for all p > 3 and any inclination. The error is larger for points on the far side of the ring (0 < ф < m, y > 2/2). 
However, even at an inclination angle of 60^, the accuracy is quite adequate, as shown by the right panel in Figure 11. 


B. TRANSFORMATIONS OF FIELD COMPONENTS 


In the analysis given in the main text, we assumed that the magnetic field components B,, By, B; are specified in the fluid 
frame. Under the usual assumptions of ideal MHD, the electric field vanishes in this frame: E, = Ey = E; = 0. Alternatively, we 
might wish to work with field components in the P-frame: po? 3 BO’, B. EO. ЕФ Я EY (the electric field does not vanish in this 
frame). 

The two frames are related by a Lorentz transformation with velocity B (expressed in terms of B and y, see eq 8). The transfor- 
mation is most transparent when we rewrite the radial and tangential field components in terms of “parallel” and “perpendicular” 
field components relative to the velocity: 


B® = cosy B® +siny BY, В --siny BP + cosy By”, MED 
B® = cosy BP -siny BY, By = siny вр + созу BY”, (В.2) 


with similar expressions for E'P and B. The transformation rules are then 


E, = ЕФ (B.3) 


_ pP) 
By = By I? 


POLARIZED IMAGE OF A RING ORBITING A BLAcK Нов 25 


B, =yBP +ByE®, В, = y BP - ByES, (B.4) 
BP-yB,-ByE, B® =yB,+ByE., (B.5) 
E, = y EP - gy EP, E, = y ED + gy B®, (B.6) 
ЕФ =yE,+ByB, ЕФ - y E, - ByB,, (B.7) 


where, as usual, y = (1 — 2)-!2, 
Using the above transformations, if we are given B,, By, В, in the fluid frame, we can solve for B®) and E® in the P-frame: 


B®) = (cos? y + y sin? y) B, – (y — 1) cos y sin y Ву, (B.8) 

BO =-(y - 1) cos y sin y B, + (sin? y + у cos? y) By, (B.9) 

BO =уВ,, (B.10) 
ЕФ) = -By sin y В., (B.11) 
ЕР =By cos y В., (B.12) 

E? = By sin y B, — By cos y Bg. (B.13) 

Similarly, if we are given the magnetic field components in the P-frame, we can solve for the other field components: 

B, = [cos? y + (1/y) sin? y]B® + ((y — 1)/y) cos y sin y BP, (B.14) 
B, - ((у — 1)/у) cos y sin y BP + [sin? y  (1/y) cos? y] Bj, (B.15) 
B, -(/y) BO), (B.16) 
ЕФ =- sin y BO, (B.17) 
ED = cos y BP, (B.18) 
E? = BsinyB? – B cos xB”. (B.19) 


These transformations are provided here for the convenience of readers who might prefer to work with field components in the 
Schwarzschild frame. 


С. EMISSION LOCATION VS. OBSERVED COORDINATES 


The radiation emitted by the point P in the ring at (R, ф) reaches the observer at sky coordinates (x, y), which we can write in 
terms of polar coordinates (p, p) as described in equation (28). Here we work out the relation between these two coordinates. 
The relation between ф and ф is straightforward. Since the observer frame is tilted with respect to the ring plane by a rotation 
angle 6, around the line of nodes, and since the geodesic lies entirely on a plane (because we have limited our analysis to the 
Schwarzschild spacetime), we find 
tan ф = tan $ cos 6s. (C.1) 


This relation can be used to translate ¢ to y and vice versa. For the analysis in Appendix D, it is useful to express q in terms of ф 
up to quadratic order. The corresponding relations are 


sing > sing + (1/2) sin? 0, sin ф cos? Q9, cosó —> cosg - (1/2) sin? 0, cos o sin? g. (C.2) 


To calculate the mapping between R and p, consider the G-frame (Fig. 1), where the geodesic lies in the xz-plane. At the 
emission point (x, y, z) = (А, 0,0), the geodesic makes an angle « with respect to the x-axis, where a is given by the Beloborodov 
approximation (4). Since the angular momentum around the y-axis in the G-frame is conserved, we have 


Кѕіп с 


(1 _ 2 


p = ky = КК = (C.3) 


Squaring both sides, 


R?(1 — cos? 
giu Abc Eee ae Gro d av ОВИ aie 48st d dad) (C4) 


@ Б 
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This directly gives p in terms of R and ¢; conversely, the quadratic equation can be solved to obtain R for a given p and ф. 
Equation (C.4) is exact, except for the fact that we used the Beloborodov approximation (4) for cos a. 
Since (0p/OR), = 0, the Jacobian determinant |J|, which describes the transformation of differential area elements between 


(К, ф) and (p, q), is given by 


др 


ІД = (2) E = 1 [к + 1) - (А – D) sin? в, sin? +2 sind sin | 
ф R 


2 
sec” Ø сов bo | (С.5) 


OR} ,\RO¢ 1 + tan? $ cos? 6, 


D. SERIES EXPANSION TO QUADRATIC ORDER 


The analysis in sec. 2 is exact, modulo the Beloborodov approximaton, and is convenient for numerical calculations. However, 
for analytical studies, we need simpler relations. For this, we expand all the equations up to second order, treating the quantities 
sin 6, B and 2/R, which describe tilt, relativistic velocity and gravity, as being small.’ The relevant series expansion results are 
given below. In each equation, the second-order terms are shown inside square brackets. 

The observed coordinates (х, y) of the geodesic emitted at location (А, ф) in the ring are given by 


1 R 
х=(К + 1) соѕф+ -5g 0050 + 2sin A sing соз ф — 7 sin? 6, зіп? g cos o А (D.1) 
р 1 А i «2 R - 2 +3 
у=(К + 1)sing + ор Sn + 2800 Gp sin 9- Sin 0,sin gl. (D.2) 


In deriving these results, we first evaluated equation (22) and then made the substitutions given in equation (C.2). The latter 
substitution is made in all the subsequent results presented in this Appendix; thus the results are expressed in terms of the 
observed azimuthal angle y. 

To quadratic order, the Doppler factor 6 is 


6=1- | (D.3) 


1 2 1 2 
RU Ё + aR T É cosy + Взїпб, sin(x + Ф) 


Note that Doppler boost due to azimuthal velocity is described by the last term, B sin 0, sin(y + y), which appears only at second 
order in the small quantities sin 0, and £. This is one of the reasons for expanding the equations to quadratic order. 
Assuming that the spectral index œ, = 1, the intensity of the linear polarized radiation at the observer is given by equation (15): 


IP] = ó* l IBP sin? z. (D.4) 


Expanding to quadratic order, the term 182 sin? С is given by 
324-2 2 : : 4 
IB sin” ё = Be, + {2 sin 6, sin(7 + ф) — R cos 7 + 2B cos(x — n) | B4B. 


"n 2 : 
+|— (si 0, sin(7 + 9) — R cos n + B cos(x — n| Be 


4. „ ; 4 a ; ; 48 2| p2 
+ -R sin б, sin o + m + sin’ 0, + 26 sin Ө, sin(y + 9) — p X +6 | В: 
ae ; 
-Е Sin 6, cos y sin PBegB-| . (D.5) 


We have written the result in terms of the parameters Ba, 7, B. of the magnetic field in the fluid frame (see eq 12). This is helpful 
for the discussion in sec. 4. Note that, in the absence of any equatorial magnetic field, the only contributions are at the second 
order (because the only terms with В? are inside the square brackets). Since the observed intensity is directly proportional to 


18р sin? Z, we need to expand to quadratic order to handle models with pure В.. 
To quadratic order, the path length /, in equation (13) is 


1 4 4 4 
=1+ 5 B+ gm + sin? 0, + 28 sin 6, sin(xy + ф) — # cosy - g sin Go sin : (р.6) 


Н 


7 Because the solution for the coordinate x involves a division by sin б, it is necessary to keep terms up to sin? 0, in the expressions leading up to this quantity. 
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We calculate the linear polarized intensity |P| as the product of the three terms, б^, lp and |В|? sin? £ (see equation D.4). This 
gives 


4 2 
|P(y)| -hı m H (В, + B;) + 2 (sin 6, cos o + B sin y) BB: + 2 |- + B cos y + sin Ө, sin e) В.В, 


Б B3 s 108 | | | p 2 
+ pesing + -z + п 6, cos 2ф + —- cosy + £ sin 6, (Sin(y — p) – 4 sin(y + ф)) — 7-(4 + cos 2y) B; 
2 6 1 6 а 
+(-2 sin б sin ф + RTI sin? 0, cos 2% — B sin 6, (4 sin(x + q) + sin(x — ф)) + T cos Y — Ea = c0s2y)} B; 
ГА 
4 4 4 
+ shinee ae eS” Г, В? 
R R2 R $ 
4 . i34 : 4B . 2. 
+ g 512 б cos ф— sin 0, sin 2y — 28 sin 0, cos(x — ф) + nx Ё sin2y | В.В, 
8 8 12 16 8 
t E sin 0 cos o — T sinx] BB; + (5 sin @ sin ф + mo + cos} вв ; (D.7) 


where we have written the answer in terms of B,, By, B. in the fluid frame. 


The electric field components Ex, Ey, which are normalized such that they are proportional to sin Z |В| (see eq 25), аге 
2 . 
Еу =- sing B, — cos p B, — | sin 0, — R sing + Bsin(y + q)| В, 


+ 


2 2 1 
s sin 6, sin? y + m sing + 2 sin? @ sin? ф + sin 0,(соѕ y — cos(y + 2¢)) 
9) 2 
-2 sin(x + y) + (sin ф + sin(2y + e» B, 
I . 1.5 В . А А 
+ -R sin @ sin 2y + 8 sin“ 6,(5 cos o — cos 3ф) + 2 sin (зїп y + sin(xy + 29)) 
В? 
4 9s y — cos(2y + e Bg 
2. ‚2 
+e sin 0, sin“ Ф В, |, (D.8) 
| 2 
E, =соѕ ọ B, — sing By + -R cos Q + B cos(x + ф) | B; 


+ 


1 2 1 
[5 sin 0, sin 24 — 2 COS Q — 8 sin? (cos ф — cos 3g) + А sin (зїп y — sin(x + 29)) 


2 
+2 cos(y + ф) — 5 (cos ф + cos(2y + e) B, 


2 1 
+ (s sin б cos? ф — 8 sin? &,(sin ф + sin 3y) — 5 sin 0,(соѕ y + cos(y + 29)) 
8° 
+ mn y — sin(2y + e) В, 
io : 
^R sin 6, sin 29 В, |. (D.9) 


From Ex, Ey, we can obtain the observed field components, Е, оь, Ey,obs, from equations (26), (27). We can then compute the 
Stokes parameters Q and U via 


О = Ea, Е? = (Е? – ED) 0 I^, U = 2E; obsEy,obs ан уа", (D.10) 


We can also calculate |P| = E s t E s but this will simply reproduce the answer given in equation (D.7). We do not write 
down the results for Q and U as the expressions are large. Instead we define the complex polarization P(y) in the usual way (see 


28 


eq. 29), and expand it in a Fourier series as described in Palumbo et al. (2020), 


Ls 
PO) = Q9) + Uy) = 7- b3 Bnei, 


m=— 00 


(D.11) 


To zeroth and linear order there are only two non-zero coefficients, 6; and 82, and to quadratic order, there are five non-zero 
coefficients, Bo — B4. The expressions for these coefficients are given below (second-order contributions are shown inside square 


brackets): 
P ; 
Bo E sin? 6, (B? + 3B; – 4B? - 2iB,B)| 
1 ; 
= l: sin? 6, (е2 — 2) B2, + sin? eB | 


В. = 251п 0, (-iB, + Bs) В, + 


= –2іѕіп бе"! ВВ, + 


j | 3, 
Е (2 - е") sin 6, + if sin 6, ^ + ser) B2 
4i | 10i | 
+(— - 2ige'* | sin 6B2 + — sin б„е'ВВ,|, 
R CUR à 
4 EV. T 2 | 
В=-|1- 5 (В, +iBy) –2| вей - = т) (Br +1В,)В. 


+ 


2 iB 
E a siny — 10icos y) + pg ; i8 + en) В? 


4 
x x pe "Js 2 +(- > + ву Ai cos y) + 418 ) В.В, 
б 


LI. " eX | В,В. + 0 UP gm: 
2 R R z 


2 
--(1- s) eg. + (же ES B. 
R 


+ 


+ 


"(m ond +9) В» 


2 | | | 
«dz - agi) + É 8 (e 2iy cs дет) – Ё (e (6 cosy + 2e*) + 2e) Be 
4 4 | А. 
"NL: et — pox) В? a [- em + 8B ач) В.В, |, 
m R z R2 R 


rm [ _ Fer sind, (в, " iB) _ Z пө, (В, + В,) a, 


For each £,,, coefficient, we give the result both in terms of B,, By, Bz, and in terms of Beg, rj, Bz. 
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